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Abstract 

This paper deals with the study of collapsing plane symmetric 
source in the presence of heat flux. For this purpose, we have calcu¬ 
lated the Einstein field equations as well as Weyl tensor components. 
The conditions for the conformal flatness have been determined. The 
interior source has been matched smoothly with the exterior geometry 
in single null coordinate. It has been found the pressure is balanced 
with the out going heat flux and the continuity of the masses in two 
regions has been noted. A simple new model of collapse has been 
proposed which satisfies flatness condition, also we have discussed the 
physical properties of the model. For our model, we have calculated 
the temperature profile by using the approximation scheme. 
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1 Introduction 


The concept of gravity in terms of spacetime curvature was introduced for 
the first time by Einstein. According to him, curvature in the spacetime is 
produced due to the presence of matter in which it exists. This curvature 
is taken as a source of gravitational held. He formulated a set of equations 
known as Einstein held equations. The solution of these equations have space- 
time singularities. These singularities are formed by gravitational collapse. 
Gravitational collapse is one of the important problem in general relativity. 

In astronomy a gravitational collapse is a phenomenon in which a massive 
star dozens of times larger than the solar mass contracts under the inhuence 
of its own gravity. It occurs when the internal nuclear fuel fails to supply 
sufficiently high pressure to counterbalance the gravity. There are two types 
of spacetime singularities, one is a black hole and the other is a naked. The 
singularity covered by an event horizon is called black hole and an uncovered 
singularity is called naked. Event horizon is the boundary of the black hole. 

Many researchers investigated the solutions of Einstein field equations 
with heat flux (Schafer and Goenner 2000, Ghan et al.2003, Herrera and 
Santos 2003). Ghosh and Deshker (2007) investigated exact non spherical 
radiating collapse using junction conditions. Ahmad et al. (2013) investi¬ 
gated the gravitational collapse with heat flux and gravitational waves. 

Initially Bronnikov and Kovalchuk (1983) studied cylindrical symmetry 
and found some exact solutions. The collapse of dust spheroid was studied 
numerically by Shapiro and Teukolsky (1991) using cylindrical spacetime. 
They concluded that the end product of collapse was either a black hole 
or a naked depending on the spheroid compactness. Barrabes et al.(1991) 
investigated the exact solution for collapsing convex shell. They concluded 
that in some cases they found the absence of apparent horizons. Sharif and 
Ahmad (2007) studied high-speed cylindrical collapse of two perfect fluids. 

Gutti et al.(2012) investigated the gravitational collapse of an inhnite 
cylindrical distribution of timelike dust, using the matching conditions. The 
part of anisotropy, radial heat flux and electric charge over the dynamics of 
non adiabatic cylindrical collapse by using the coupled equation was studied 
by Sharif and Abbas (2011). Abbas and his collaborators (abbas 2014a, 
Abbas 2014b Mahmood et al.2015) investigated the expanding and collapsing 
solutions for charged plane and cylindrical stars. The shearfree cylindrical 
collapse was studied by Di Brisco et al.(2009). Herrera et al. (2004) studied 
spherical shearfree radiating collapse and conformal flatness they developed 
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a relation between dissipation and density inhomogeneity. Also, Herrera and 
Santos (1997) explored the properties of anisotropic self-gravitating spheres 
and discussed their stability using the perturbation method. Herrera and his 
collaborators (Herrera et al.2008a, Herrera et ah 2008b, Herrera at al. 1989, 
Herrera et al. 2009a, Herrera et al. 2010, Herrera et al. 2012) have discussed 
the stability and applications of anisotropic solutions to stellar collapse. 

In this paper, the work done by Herrera et al. (2004) is extended by taking 
the plane symmetric gravitational with heat flux. The plan of the paper is 
as follows: The matter source and held equations are presented in section 2. 
In section 3, we provide the formulation of solutions by taking into account 
conformal hatness. In section 4, we construct a simple dissipative model 
and studied the temperature prohle of the proposed model. We conclude the 
results in the last section. 


2 Interior Matter Distribution and Field Equa¬ 
tions 

Assume that a given four dimensional plane symmetric spacetime is divided 
by a three dimensional hypersuface S into two regions interior and exterior 
denoted by H“and respectively. We consider a particular plane symmet¬ 
ric spacetime in interior region given by 

ds't = z)dt^ + z) [dx^ -|- dy'^ + dz"^]. (1) 

Here the coordinates are labeled as follows: = t,x^ = x,x‘^ = y and 

x^ = z. We assume that the energy momentum tensor in the interior region 
is given by 

Taf} = id-+ P) WaUly + Pgay + + QaWiS, (2) 

where fi is the energy density, P the pressure,re" the four velocity of huid 
and g" the heat hux satisfying gQ,rc" = 0. For interior spacetime ([I]) in 
comoving coordinates, we have 

(3) 

9“ = (4) 


where g depend on t and z. 
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In this case the rate the rate of expansion scalar 0 = is given by 


0 = 3 


AB' 


Weyl tensor CajS'ys has the following non-zero components for metric 
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For the interior metric ([1]) the non vanishing components of Einstein held 
eqnations are given by 



GS, = -2 (ll) / =(11) 


Using Eqs. llTT]) and ([5]), we get 

kqB'^ = ^ 0 ,^, ( 12 ) 

This indicates that for g > 0 heat hux is directed outward, implying that 
0^2 > 0. For g = 0, Eq. (IT^ implies that 0 ^ = 0 which shows that collapse is 
homogeneous (Herrera and Santos 2003). The Taub’s mass function (Zannias 
1990) for the plane symmetric spacetime is given by 


m{t, z) 


(gii)^ 

2 


pl2 _ 
“12 ~ 


B 


A^ 



(13) 


4 



















Differentiation of m(t, z) with respect to z and t along with held equations 
dSD-dllD gives 


m {t, z) 
rn{t, z) = 


PB^B^t + 


(14) 

(15) 


The above equations show that the heat hux may effect the gradient and 
time derivative of m(z, t), and the dissipation reduces the amount of matter, 
hence the rate of collapse slows down probably. 

The scalar of the Weyl tensor in terms of Kretchsmann scalar = 
Raises, the Ricci tensor and Ricci scalar R is dehned as 

= R - 2R^^Rc,f} + ^R^. (16) 

o 


Using Eqs.(l6]), (IH]) and flTdl) in Eq. flTHl) . we obtain Weyl tensor in terms of 
pure gravitational mass given by 





KjJ, 

T. 



(17) 


where m,, = m — 


3 Conformally flat solution 


In this section, we impose conformal hatness condition to metric ([T]), this 
requirement implies that all the components of the Weyl tensor must be 
equal to zero. Hence, from Eqs.([n]) and ([7]), we see that if C 1212 = 0, this 
condition is justihed and we have 

^22 Bzz\ B ^ BA 

hr-ir)-2if(T-if) = “- 

Integrating above equation and after reparamatrizing t, we get 



A=[C^{t)z + l]B. 


(19) 


where Ci is arbitrary function of t. Also, from the isotropy of pressure 
Gii = G 33 with Eq. lfT^ . gives 


B 


0 , 


( 20 ) 
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integration of Eq. fl2U]) yields 


( 21 ) 


C 2 (t)z + C 3 {t)' 

where C 2 and C 3 are the constant of integration depends on t. Stephani 
(1967) and Kramer et ah (1981) investigated all conformally flat solutions 
with g = 0. 

Conformal flatness imposes C = 0 and put ttt-c = 0 in Eq. flTTl) . we get 

m = (22) 
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Using Eqs. flTl)) and fl2^ we get 

/i' = gE^e, (23) 


which describes that for g > 0 and 0 < 0 then p' < 0 implying that the 
density reduces with increasing z, while from (E]) and flT^ we get 



v' = ^(0%, 

(24) 

integrating (|2^ we get 

02 

(25) 


Kg = — + g{t), 

where g depend on t only. 




Putting solution (IT^ and (|^ into (IHl), (E]) and (ITTil it follows 
- - 

KzP = —-— [2((72Z + C‘^){C2Z + C 3 ) — ‘i{C2Z + 

LjiZ + 1 

2 C 

— —{C 2 Z + C‘Y){C2Z + C‘Y)z\ 

OiZ + 1 

+TT^[CiCIz + ZCl - 2C1C2C3], ( 27 ) 

Kg = 2(CiCs-C’2)(^g^^)7 (28) 
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The expansion of the flnid given by ([S]) along with Eqs. lfT^ . fl?T]) and (125]) . 
yields 


C 2 Z + Ca ^ • nqz Ciz + 1 

Ciz+l ^ I “ 2 (Czz + Czf 


(29) 


The above eqnation show that for g = 0 the contraction is homogeneons, 
however for q ^ 0, collapse is inhomogeneons, which has already been con- 
clnded in ([12]). From Eqs. fl26|) and (|25|), we see that g{t) = —4(71. Here, we 
matched the interior spacetime with the exterior spacetime that is described 
by plane symmetric spacetime in single nnll coordinate given by 


ds‘^ = — 2dudz + z^{dy^ + dz^)^ (30) 

z 

from Eqs. ([ID and (I5UD on hypersnrface S along with field eqnations ([S P - lfTTl) 
and the mass fnnction ffTSp . we get 
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Using Eqs. (I27D - (128D and fl3^ we obtain 
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4 A simple model 

A simple approximation solution satisfying the junction condition fl34l) for 
the functions (7i(t), (72(t) and Csit) is 

Ci = 6Ci(t), C2it)=0, C3{t) = ^, (35) 
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where 0 <e<<l, 0 <t<—cxo and a > 0 a constant proportional to the 
total mass in the plane. Using fl5BD and fIMD we get npto order of e, 


Cl + — Ri 0, 


integrating fl36p we get, 


Ci{t) Ri Ci(0)exp[-]. 

■Zs 


Pntting the solntion fl35|l into fl26|) - fl28|) . we get 

12 a 2 
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above conditions hold for a realistic process in stellar evalnation. We wonld 
like to mention that for e 7 ^ 0 the valne of t is constrained by some physical 
sitnation. Thns, if one reqnires that the central pressnre remain less than 
the valne of the central energy density, then there exists the ineqnality 

^ > eci. (41) 

t2 

From Eq. fl55]l and Eq. flS^ . we see that the energy density and pressnre 
decrease in the onter region dne to dissipation. It can be noted from (I40|) . 
that there is increase in the heat flow ontwardly. 

From flT^ with (IT^ and fl2l]) we get 

2 

m{t, z) ^-{1 — 2eciz), (42) 

a 

which indicates that mass of the gravitating system is being decreased dne 
to heat flnx. Now the expansion scalar fl35l) becomes 


0 Ri ^(1 - eciz). 
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(43) 





The above equation shows that dissipation causes to decrease the rate of 
collapse. The adiabatic index in the presence of heat flux is defined as follows 


/p\ 

dln^ yPy \A/ ' 


(44) 


It gives the dynamical instability at a particular time, now computing (144)) 
at 2 : = 0 and z = z-s (using Eq. (l35|) and (|37l) ) upto order of e, we get 
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The above equations imply that the collapse is homogeneous at the center 
and on the hyper surface S. 


4.1 Calculation of Temperature 

This section deals with the calculation of temperature prohle T{z,t) for our 
model. The heat transport equation in Maxwell-Cattaneo theory (Herrera 
et al.2004) is 

+ q^ = -Kh^^iTp + apT), (47) 

where the quantities r, K and h°‘^ = g°‘^ + w°‘w^ are the time relaxation, 
thermal conductivity and the projection tensor respectively. Considering 
Eqs. ([I])-Oil) then (|471) becomes 

TB{Bq)^t + qAB^ = -K{TA)^,, (48) 

Applying Eqs. ([19]), (l 2 T|) and (l28l) into (1481) and taking C 2 = 0, we obtain 
upto 0 (e) 

kK 

T{CjCsC3),t + O 1 O 3 = [T(ci^+ 1)],,. (49) 

Now if there is no dissipation we have (Oi = eci = 0) from (149|) it follows 
that T = To{t) this shows that for non-dissipative systems, the temperature 
is homogenous through out the system. Thus for dissipative case Ci ^ 0, we 
have the following form of temperature prohle 

T = TQ{t) + eT^{z,t). (50) 
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substituting flSU|) into fl^ we get upto 0(e) 


^ , /4a 



4recia^ 

nKt^ 



(51) 


Here, it is assumed that X > 0 and Tc is the central temperature. The 
dissipation can decrease the temperature at the loss of energy, as it is evident 
from second term of term of above equation while last term predicts the 
support of relaxational effects. 


5 Conclusion 

It is well known that gravitational collapse of a stellar object is extremely 
dissipative process. Therefore, it becomes important to study the effects of 
dissipation during the collapse of a radiating star. During the dynamical 
process of gravitational collapse of a non-adiabatic star, the enough amount 
of energy losses in the form of outward heat flow and radiations. Herrera et 
al. (2009a) derived the dynamical equations by including the dissipation term 
in the source in the form of heat flow, radiation, shear and bulk viscosity 
and then coupled the held equations with the causal transport equation. 
The inertia due to heat hux and its effect during the dynamics of dissipative 
collapse with outgoing radiations was studied by Herrera (2006). 

In this paper, we investigate the plane gravitational collapse of radiating 
huid with conformal hatness condition. For this purpose, we have explored 
some the effects of dissipation on the dynamics of self-gravitating body. By 
determining the conformal hatness condition of the collapsing star, we have 
found that the system of equations is exactly solvable and there appears 
three arbitrary integrating functions which depend on t. This implies that 
conformal hatness conditions to leads to the homogeneous solutions. Further, 
the interior solution has been matched smoothly with the plane symmetric 
metric which is dehned in single null coordinate. We would like to mention 
that exterior geometry emits the radiations in the form of outgoing radiations. 
We constructed a simple model which satishes the matching conditions. 

We have calculated the adiabatic index which shows that gravitating sys¬ 
tem is instable at particular instant of time. The perturbation analysis of 
the adiabatic index implies that the collapse is homogeneous at the center as 
well as on the hypersurface S. Using the Maxwell-Cattaneo heat transport 
equation, we have calculated the temperature prohle of the radiating star. 
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The importance of the proposed model lies in the fact that it presents the 
influence of relaxational effects on the temperature prohle in a sophisticated 
way and thereby on the evolution of the system. It is concluded that the in¬ 
homogeneities in the energy density are directly related to dissipation, even 
though the underlying spacetime is conformally flat. 
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